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Abstract 

a 

The nonlocal Polyakov-loop-extendcd Nambu Jona-Lasinio (PNJL) model is further im- 
" proved by including momentum-dependent wave-function renormalization in the quark quasi- 

particle propagator. Both two- and three-flavor versions of this improved PNJL model are 
discussed, the latter with inclusion of the (nonlocal) 't Hooft-Kobayashi-Maskawa clctcrmi- 
nant interaction in order to account for the axial U(l) anomaly. Thermodynamics and phases 
4^ are investigated and compared with recent lattice-QCD results. 

Oh 

§^ 1 Introduction 

The Polyakov-loop-extended Nambu and Jona-Lasinio (PNJL) model has become a widely used 
"minimal" approach to deal with important aspects of the QCD phase diagram, in close contact 
and comparison with lattice-QCD results whenever available. Its basis is the NJL model [1-5], 
— . generalized by incorporating Polyakov-loop dynamics that is controlled by a suitably chosen 

VO effective potential [6-10]. Once the dynamical input is fixed in terms of a few parameters 

^ by reproducing basic facts from hadron physics in vacuum and QCD thermodynamics, the 

^ PNJL model permits to explore, at least schematically, regions of the phase diagram that are 

not accessible to lattice QCD but nonetheless of principal interest. A particularly attractive 
y—i feature of this model is its ability, through the coupling of the Polyakov loop to the quarks 

^ as dynamical quasiparticles, to study the correlated pattern of the deconfinement and chiral 

crossover transitions. 

In recent work [11-16] the PNJL model has been generalized and extended by introducing 
nonlocal effective interactions of the quark fields. This was an important step in order to avoid 
the traditional sharp momentum-space cutoff characteristic of former (local) NJL models. The 
nonlocal version now permits to establish contacts with QCD at the level of quark quasiparticle 
properties. In Refs. [11, 12] the first part of such a program has been accomplished by gener- 
ating momentum-dependent dynamical quark masses from a generalized gap equation, in close 
correspondence with (Landau gauge) Dyson- Schwinger equations or (extrapolated) lattice-QCD 
results. At this stage, however, quark wave-function renormalization effects were not yet incor- 
porated. The residue of the quark propagator at its quasiparticle pole was strictly Z = 1. The 
commonly used chiral quark interactions are not explicitly momentum dependent and, conse- 
quently, do not provide a mechanism for reducing the wave- function renormalization factor Z{p) 
from unity (at least not at mean-field level). 
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In the present work we generalize the PNJL model further in order to improve on this issue. 
A vector-type derivative coupling between quarks will be introduced such as to reproduce results 
for Z(p) from lattice QCD. The two-flavor version will be dealt with in Sect. 2, followed by the 
three-flavor case with inclusion of the axial U(l) anomaly in Sect. 3. The thermodynamics 
based on this extended PNJL model will be derived in Sect. 4. The discussion of chiral and 
Polyakov-loop crossover transitions at zero quark chemical potential will be updated, also in 
view of recent lattice results. Finally, extrapolations to nonzero quark chemical potentials and 
the resulting phase diagram will be explored. We present a summary and an outlook in Sect. 5. 



2 Generalized nonlocal Nambu— Jona-Lasinio model 

2.1 Quark propagator and wave- function renormalization 

Consider as a starting point the full, renormalizcd Euclidean quark propagator S(p), the inverse 
of which is represented by the Dyson-Schwinger (DS) equation 



+ 




S{p)-^ = -Z2{p-mo)-Z, I ^,g^DZ{p-q)^l^S{q)^r,{q,p). (2.1) 

It involves the renormalizcd (dressed) gluon propagator D'^ and the full quark-gluon vertex func- 
tion 11/ , together with renormalization constants Zi,Z2. The solution of the DS equation (2.1) 
has the general form 

The functions A{p'^) and B{p^) depend on the renormalization scale and encode the Lorentz- 
vector and -scalar dressings of the propagating quark by its interaction with the full gluon 
field. The wave-function renormalization factor, Z{p^) = A~^{p'^), accounts for the deviation 
of the quark's quasiparticle strength from unity, induced by momentum-dependent vector-type 
interactions. The dynamical quark mass function, Af(p^) = B{p'^)Z{p^)^ is generated primarily 
by nonperturbative mechanisms involving the scalar quark condensate (V'V') that acts as an 
order parameter of spontaneously broken chiral symmetry. In the chiral limit (mo = 0), the 
mass M(p^) remains strictly zero at all orders of QCD perturbation theory. 



2.2 Nonlocal Nambu— Jona-Lasinio model with two quark flavors 

In Refs. [11, 12] we have already shown how a nonlocal generalization of the Nambu and Jona- 
Lasinio (NJL) model [1] can be constructed starting from the SD equation (2.1) in the case 
oi Z{p^) = 1. We are now going to write down the extended nonlocal action, 5e, for the 
two- flavor NJL model including a momentum-dependent renormalization function Z{p^) (fol- 
lowing Refs. [14,15,17-19]). Assuming separability as in Refs. [11,12], i. c, g'^D^{p~ q) — )• 
^ C{p) C{q) Sjj^i^Sab, where G is a coupling strength of mass dimension —2, we can write in Eu- 
clidean space-time: 

Se= I d^x U{x) [-10 + mo) ^l:{x) - ^ \ja{x)ja{x) + J{x)J{x)] I , (2.3) 
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with the quark fields ip{x). We work in the isospin Hmit with equal current quark masses for up 
and down quarks, = = mo- The nonlocal currents ja{x), J(x) are given by 



a{x) = j dSC(z)VS(x + r„V'(a;-|) , (2.4) 

J{x) = |dS^(z)V;(x + 0^v(x-|) • (2.5) 

Here the usual chiral scalar-isoscalar and pseudoscalar-isovector NJL operators, /lie{o,...,3} = 

_ -f> 

(1,175-?) are now supplemented by vector currents with derivative couplings, il){x')$ip{x) := 
V'(.x') {$ijj){x) — {$'i[j){x') ip{x). The nonlocality distribution C{z) is the same as in Refs. [11, 12]. 
In momentum space: 

'exp [-p'^d?(j/2] for < 



const.. ^ forp2>A^ ^'"'^ 



p2 



normalized as C{p = 0) = 1. The high-momentum behavior, C{p'^) oc as(p^)/p^, with the QCD 
running coupling Q!s(p^), is known from the operator product expansion for the dynamically 
generated quark mass [20] . It is matched at a scale A to the nonperturbative infrared sector of 
C(p^) that is conveniently parametrized as a Gaussian. Its width reflects a typical instanton size 
or the characteristic gluonic field strength correlation length, dc ~ (0.3-0.4) fm. The matching 
scale A is fixed by requiring continuity and differentiability of C{p^). The detailed low- momentum 
form of (2.6) is actually not essential because all relevant integrals involve C(p^) multiplied by 
powers of p^. 

All four quark currents jai^) with a = 0,1,2,3 have the same distribution C{z) to ensure 
chiral invariance. The new element in this generalized nonlocal approach is the vector-type 
derivative current J{x) of Eq. (2.5) that is separately invariant under global chiral transforma- 
tions. We follow Refs. [18, 19] at this point. The derivative coupling proportional to Judf^ in J{x) 
induces the momentum-dependent wave-function rcnormalization in Eq. (2.2). Apart from the 
overall coupling G that has been factored out in Eq. (2.3), the coupling strength now features an 
additional independent parameter k with mass dimension 1 in the denominator at each vertex. 
The nonlocality distribution T{z) associated with J{x) is again normalized in momentum space 
as T{p = 0) = 1. It is given a Gaussian parametrization, 

.F(p2) = exp [ V4/2] , (2.7) 

with dp to be adjusted such that the momentum-dependent renormalization factor Z(p) is well 
reproduced. Note that this ansatz is actually consistent with the one for C{p) in Eq. (2.6): 
according to Ref. [21] there is no modification of Z{p^) = 1 for p^ > A^ coming from a Landau- 
gauge operator product expansion up to first order in the current quark mass. 

Given the action (2.3), we apply the standard bosonization procedure. Consider the partition 
function 

Z = j (2.8) 

and introduce five bosonic fields (j,v and vr^ plus five auxiliary fields S,V.Pk {k G {1,2,3}). 
Next, insert a ""one'" in terms of (Euclidean) delta functions on the right-hand side of Eq. (2.8), 



J 9S m 9Pk 5{S - jo) SiV - J) S{Pk - jk) 



= I &Pk &v ^TTfe exp <^ / d^z 



{/■ 



(^{S - jo) + V{V - J) + TTkiPk - jk) 
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and write the partition function as 



exp I - y d^rc [ip {-i^^'d^ + mo) + o-jo + vJ + -KkJk] | (2.9) 
X exp < / d^x 



G . 



{(rS + vV + TTkPk) + ^{S^ + V^ + PkPk) 

The inserted delta functions imply a replacement of jo, J and jk by S, V and P^, respectively. 
The first exponential, after re-insertion of Eqs. (2.4) and (2.5) for the currents, has terms of the 
form ■0^^. The path integration over the fermionic fields and ■0 is carried out in the standard 
way: 



/ 



Z= \ 9aSsv ^TTfe det 



J &S m 9Pk exp I J d'^x 



{aS + vV + TTkPk) + ^{S^ + V^ + PkPk) 



(2.10) 



The integration over the auxiliary fields S, V, Pk is readily performed by completing the square, 
leading to 

Z=J\f J ^a^v^TTk det^/exp|-^yd'^x[a^ + ?;^ + 7r2]| , (2.11) 

where det is the fermion determinant. In momentum space one finds after Fourier transfor- 
mation^ 



p + p'\ p + p' 



P-^^ ) [cr{p -p')+i 75Tfc7rfc(p - p')] 



2k 



-v{p-p') 



(2.12) 

Finally, comparison of Eq. (2.11) with the definition (2.8) of the partition function gives the 
bosonized Euclidean action, 

1 f d^ 

5bos = - lndet^(p,p') + ^ y H-pXp) + v{-p)v{p) + 7f(-p) • 7f(p)] , (2.13) 

where we have used integration in momentum space^ (for later convenience). 
2.3 Taylor expansion of the action 

In order to derive formulas for physical observables, such as the pion mass and the pion decay 
constant, a Taylor expansion of the bosonized action (2.13) is performed around the expectation 
values of the fields a, f, vf. We assume that a and v have nontrivial translationally invariant 
mean-field values a, respectively, while the mean-field values of the pseudoscalar fields tt^ 
{k = 1, 2, 3) vanish. In coordinate space: ct{x) = a + 6a{x),v{x) = v + dv{x) and 7t{x) = 6Tr{x). 
This permits to write the action (2.13) as 



(2.14) 



^We use C{p) = / d*zC{z) e ^^'^ , <j){p) — f d'^z cj>{z) e ^'''^ , where cj>{z) stands for an arbitrary field G {a, v, if}. 
We omit the "hats" on the Fourier transforms and set henceforth C{p) = C{p),d>{p) = (f>{p} . civ. 

^The arguments of the quadratic expressions in a, v, n follow immediately from the Plancherel theorem under 
the assumption that those fields are real-valued in coordinate space. 
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The mean-field approximation is defined by the requirement of a vanishing first derivative of the 
action »Sbos) 

=0, ior ^e{a,v,Tr}. (2.15) 

{(7=d-,v=v,n=0} 



dip 



The mean-field action, , is thus obtained by replacing (j{x) — )• a,v{x) — )• v,tt{x) — > 0. In 
momentum space this implies the substitutions a(jp) — >■ (27r)^(5(^)(p), v(jp) — >■ v {2itY5^^^(jp) 
and, of course, t({p) 0. With this and using^ Indet = Trln, we find 



d^y 
(27r)4 



Trln [Z-^{p){-f + M{p))\ + 



-I- 
2G 



(2.16) 



where F'-^-' denotes the four-dimensional Euclidean space-time volume and the number of 
colors (iVc = 3). We have introduced the wave-function renormalization factor, Z{p), 



and the dynamically generated constituent quark mass, M{p), 

M{p) = Z{p) {mo + aC{p)) . 



(2.17) 



(2.18) 



The quadratic term 5^^^ follows from the second-order term of the Taylor-series expansion 
in 5(j): 



1 



d^jo 
(2^ 



[Ga{p^) Sa{p) 6a{-p) + g^{p^) 6v{p) Sv{-p) 



(2.19) 



+2g^,v{p^) Sa{p) 5v{-p) + g^ip") S7r{p) 67r{-p) 
The inverse propagators for the a and v fields are given by^ 



5a{p) da{—p) 
= iz-8Nr, 



G 



{cr,i',7ri=0} 



e\q)Z{q+)Z{q 



M{q+)M{q- 



{q+f + M^{q+) {q-f + M^{q- 



gvip" 



ga,v{jp^ 



bos 



dv{p) Sv{—p) 



1 8Nc 



{a,v,0} 

^^q':F\q)Z{q^)Z{q- 



q+-q-- M{q+)M{q-) + 




'f {q-f - (g+ • <l-f 


/(V) 




\q+f + M^iq+) 




\q-f + M^q-) 





da{p) Sv{—p) 
n J (27r) 



{a,v,0} 

-X{q)F{q)Z{q+)Z{q- 



q-[q-M{q+) + q+M{q-)] 



{q+f + M^{q+) {q-f+M\q 



The functions In, dot and Tr extend over Dirac and flavor space only; operation on color space has already 
been applied, leading to the factor A*'c. 

"^Notice, that from the definition of the gamma matrices in Euclidean space-time, {7/a,7i/} = —26^v, one has 
tr{7^7^} = -4(5^1. and tr{7^7^7p7a} = 4 {5^v5p„ - 5^p5^„ + 5f,a6„p). 
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with = The fields a and v obviously mix. In the present context this mixing is of no 

relevance and we will not discuss it further. More important in the present work is the inverse 
pion propagator, Q-k{p'^)^ that is used for the determination of the pion mass: 



1 f d^q 



{a,v,-!ri=0} 



G 



8N, 



C\q)Z{q+)Z{q-) 



q~ + M{q+)M{q- 



(2.20) 



{q+f + M^{q+) {q-f + M^{q-) 



Note that no summation is implied in this second derivative: with tr{rj ■ t^} = 26ik only the 
diagonal elements in isospin space survive. 

2.4 Mean-field approximation and chiral condensates 

As mentioned earlier (see Eq. (2.15)), the mean-field values a and v are determined by imposing 
a minimal action 5^^. From Eq. (2.16) we obtain the following two gap equations: 



a = 8NcG 



8NcG 



(2vr) 



dV 



p2 + M2(p) ' 
p2 Z{p) 



K .1 (27r)4-^^^^p2 + j^2(p) 



(2.21a) 
(2.21b) 



with M{p) and Z{p) determined self-consistently by Eqs. (2.17) and (2.18). 

The chiral condensates can be calculated using the Feynman-Hellmann theorem by calculat- 
ing the variation of = exp [— ^S''^^] with respect to the corresponding current quark masses. 
The definition of the chiral condensate (qq) = (uu) = (dd) is^ 



(qq) = -itr lim SF{x,y) - S^\x,y) 



(2.22) 



where Sp^x,!/) is the full quark propagator while Sp \x,y) denotes its perturbative part. For 
nonvanishing current quark mass, uiq ^ 0, the perturbative contribution needs to be subtracted, 
leading to 

Z{p)M{p) rriQ 



(2.23) 



4- jvf2 [p) p^ + mQ^ 
2.5 Pion niciss and pion decay constant 

Prom the inverse pion propagator Q^nip^) the pion mass is determined by the condition 

G^{-ml) = ^. (2.24) 

Furthermore, the square of the quark-pion coupling constant, g^qq^ is defined as the residue of 
the pseudoscalar-isovector quark-antiquark amplitude, 



QQ 



p^ + 



hence 



-2 



dp2 



(2.25) 



(2.26) 



^Here tr denotes the trace over Dirac space only. Since the condensates of each flavor are considered separately, 
the trace over flavor space must not be taken. 
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Finally, comparing the definition of the (inverse) pion propagator, Eq. (2.20), with Eq. (2.25) 
the renormalized one-pion state, TTk, is defined as 

TTfe = 9nqq'^k , (2-27) 

such that 5^<Sbos/(<^7r(p) STr{—p)) = + ml.. 

The derivation of the pion decay constant in the presence of the wave-function renormal- 
ization factor, Z{p), is considerably more involved than the simpler case discussed in Refs. [11, 
12]. The decay constant /tt is defined as the matrix element of the axial current Jx^{x) = 

'ip{x)^^^^^'ilj{x) between the vacuum and the renormalized one-pion state at the pion pole, 

(0|Jl,^(0)|*fe(p)) = i<5,fep^/.. (2.28) 

In order to calculate this matrix element one has to gauge the nonlocal action (2.13) by in- 
troducing a set of axial gauge fields A^^{x),i G {1,2,3}. Since we are dealing with a nonlocal 
theory this gauging requires a generalization of the parallel transport. The "connection" that 
links operators at two different space-time points x, y is given by a Wilson line, W{x, y), defined 
as 

W{x, y) = P exp I ^ ds^ 75r • A^{s) | , (2.29) 

where s runs over an arbitrary path connecting x with y. This procedure implies the replace- 
ments^ 

^(^ _ ^ W{x, X - z/2) i;{x - z/2) , 

ip\x + z/2) ^4}^{x + z/2)W{x + z/2,x). 

In the case of local operators it gives the well-known minimal gauge-coupling rule dfj, ^ + 
575r •/p (.x). 

Starting from Eq. (2.9) the replacements (2.30) enter first directly in the kinetic term and, 
secondly, in the currents jo, jk and J of Eqs. (2.4) and (2.5). The fermion determinant of the 
gauged theory reads, in coordinate space, ^ 



{x,y) 



^ (-1^2/ + ^75r • i + mo^ 5^^) {x - y) 
+ C(. - ,),o W fx, ^?l±ll ,0 a f ^ W f X 



2 

+ C(x - rtioW^a:, 5^) i7075?- *(54^) w(^, 

_ - y) [K - ^.^-A.^ .0 w(.. ^) ^ .(^) >v(^,.) . 

(2.31) 

Choosing a straight line connecting the points x and y in the Wilson line >V(x, y), we have 

W{x, y)=rexp!^^J^ da 75f • /^(x + (y - x)a) {y^ - x^) | , (2.32) 

and its derivative with respect to the gauge field A]^{t), 
5W{x,y) 



A=o 2 

i=0 



/ da 75Ti 5^^^ {x + {y- x)a) {y„ - x^,) . (2.33) 
Jo 



^Note, that this replacement apphes to tp^ and not to V". This distinction is relevant in our case because of the 
75 in the Wilson line (2.29). 

'^The 7o matrices appearing in ^'^(x, y) are the artifacts of the replacement rule (2.30) for ip'^ . Since the path 
integration extends over i/j, an additional 70 has to be included because of V'^ = V'Ho7o = V'7o- 
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The matrix element (2.28) becomes^: 
/ni ji ic^w I w det^*^ 



,4=0 

t=0 



= 8iiVe(5ifcTr|^ daq^ 

+ 8iNc6ikax 

dC{q) ^/ p\ Z{q+)Z{q-) [q+ ■ q- +M{q+)M{q-)] 



2 ^ [(g+)^ + M2(g+)][(g-)VM2(g-)] 



+ 8i N^SikX 



f)Z(g+)Z(g-)M(g- 



[(g+)' + M2(g+)] [(g-)'+M2(g, 



with 

Qa = Q + - a) Qa = Q - 2^'^ + a) 

+ .P _ P 

Note that expressions containing ^J^{p) have been replaced by (cf. Eq. (2.17)) 

vdT{p) dZ-i(p) 



(2.34) 
(2.35) 



—jr{p) = Z-\p)-l, 



K dp2 dp2 



Furthermore, the mass function M(p) is given by Eq. (2.18) and includes the wave-function 
renormalization factor Z(p). 

With Eq. (2.34) the pion decay constant is now given as 

U = ipMJi,M^^(P))^-^ (2-36) 

evaluated at the pion pole = —rn^. The pion- field component with index i is singled out in 
this expression and no summation is implied. 

Finally Taylor-series expansions of the inverse pion propagator, Eq. (2.20), and of the pion 
decay constant, Eq. (2.36), recover the Goldberger-Treiman relation, 

U9^qq = ^ + 0{ml), (2.37) 

and the Gell-Mann-Oakes-Renner relation, 

mlfl = -2mo{qq) + 0{ml) , (2.38) 

at leading order in the quark mass mo. This means that the nonlocal NJL model including wave- 
function renormalization preserves chiral low-energy theorems and current-algebra relations. 



®Tr denotes the functional trace over momentum space. According to Ref. [18] tlie integrals over a can actually 
be solved analytically. 
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mo 


G 


K 


3.0 MeV 


2.37 fm^ 


5.21 GeV 



Table 1: Current quark mass uiq and coupling-strength parameters G, k. 





a 


V 






M(0) 


-0.288 


0.493 


-1.928 


0.084 


0.142 


0.362 



Table 2: Output resulting from the parameters given in Table 1. All numbers are given in GeV. 



2.6 Pcirameter fixing 

The two-flavor case discussed in this section involves primarily the coupling-strength parameters 
G, K and the (current) quark mass mo- These parameters are chosen such that the pion mass, 
~ 0.14 GcV, and the pion decay constant, ~ 0.09 GeV, arc reproduced. At the same 
time, the normalized nonlocality distributions C{p) and J^{p), associated with the chiral scalar- 
pseudoscalar and vector-derivative interactions, are constrained by the available information 
from lattice QCD or Dyson-Schwinger-equation results for the dynamical quark mass M{p) and 
the renormalization factor Z{p) (sec Figs. 1, 2). The lattice results taken here for orientation 
are Landau- gauge simulations with two dynamical quark flavors for M{p) extrapolated to the 
chiral limit [22, 23], and with 2-1-1 dynamical quarks for Z{p) (with renormalization scale set at 
3 GeV) [24]. 

With these constraints the nonlocality length scales dc and dp oi the Gaussian distributions 
in Eqs. (2.6) and (2.7) become 

dc = 0.40 fm , dF = 0.27 fm , (2.39) 

well within the expected range of gluonic correlation lengths. The matching scale at which C{p) 
turns over to its high-momentum expression proportional to as{p'^)/p^ is then A = 0.83 GeV. 

The input parameters are summarized in Table 1. The resulting output is collected in 
Tabic 2. Note that the small value of the current quark mass, mo — 3 MeV, together with 
the large magnitude of the chiral condensate {qq), is consistent with the renormalization scale 
/i ~ 3 GeV characteristic of the lattice-QCD results for M{p) and Z{p) that have been used to 
constrain the quark quasiparticlc propagator. The scale- independent product mo{qq) enters in 
the Gell-Mann-Oakes-Renner relation (2.38). 

The parameter k, associated with the vector-type derivative coupling, comes in the combina- 
tion s/GP/k where p is a typical momentum (p ^ 0.5 GeV) appearing in the relevant amplitudes 
and integrals. This coupling strength is effectively about one order of magnitude smaller than 
VG itself, indicating a convergent hierarchy in a derivative expansion of the quark quasiparticle 
interactions. 

The dynamical quark mass M{p) in Fig. 2, extrapolated towards the low-momentum region, 
is subject to uncertainties concerning its limit M{p = 0), commonly interpreted as a constituent 
quark mass. Our standard choice is the solid line in Fig. 2 with M(0) ~ 0.36 GeV. The dashed 
curve in Fig. 2 shows an option reproducing M{p = 0) ~ 0.28 GeV (we recall that values around 
M(0) ~ 0.3 GeV are frequently used in phenomenological quark models). Such an M(0) can 
be reached with a slightly reduced coupling strength, G = 2.01 fm^, but at the expense of 
simultaneously reducing the pion decay constant to ~ 74 MeV. 
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p [GeV] p [GeV] 

Figure 1: a) The distribution C{p) (Eq. (2.6)) and b) the wave-function renormahzation factor 
Z{p) (Eq. (2.17)). Lattice data of Z{p) are taken from Ref. [24] for different current quark 
masses as indicated. 




Figure 2: Momentum dependence of the dynamically generated quark mass M{p) (see 
Eq. (2.18)). Diamonds show unquenched lattice results for the two- flavor case, extrapolated 
to the chiral limit [22]. The open triangles and circles correspond to extrapolated two-flavor 
lattice results [23] for quenched and unquenched simulations, as indicated. The solid curve is 
calculated using the input of Table 1 and C{p), Z{p) of Fig. 1. The dashed curve shows the 
effect of lowering the quark quasiparticle coupling to G = 2.01 fm^. 
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3 Nonlocal NJL model with three quark flavors 

3.1 Three-flavor action 

The action for the improved nonlocal three-flavor model is determined by the action given in 
Ref. [12], plus the extra term introducing the vector derivative interaction, Eq. (2.5): 

5e = / d^x {^{x) i-iji^d^ + m,] V(x) - f [£ix)£ix) +j^{x)j^{x) + J(x)J(x)] + 

(3.1) 



H 



where rhq is the Nf = 3 quark mass matrix, rhq = diag(mu, nid, rus). The currents are 

if(x) = j d''zC{z)i^{x+'-) Ki:(x-^^ , (3.2) 
j^{x) = j dSC(z)V5(a: + |) i75AaV(^-|) , (3-3) 

with U(l) flavor matrices Aq, (a = 0, 1, ... , 8), and 

J(x) = j d'zTiz) i^(x+'-) - I) • (3-4) 

The action (3.1) is invariant under chiral SU(3)l x SU(3)r in the limit rhq = 0. The functions 

_ ^ _ 

C{p) and J-{p) are given as in the two-flavor case; again we have 'tp(x')0ip{x) := t/j^x') {$^Jj){x) — 
{0'ip){x') ^{x). The term in the second line of Eq. (3.1) derives from the 't Hooft-Kobayashi- 
Maskawa determinant that generates the anomalous breaking of the axial U(1)a symmetry. 
This term produces an rf meson mass about twice as large as the mass of the rj meson, once the 
strength H of this determinant interaction is properly chosen. 

The bosonization procedure of the action has already been outlined in the context of the 
two-flavor model (Sect. 2.2), so we restrict ourselves to the result for the partition function: 



J ^TTa det J J 98^ 9Pa 9V exp I j d^X {aaSc + TTaPa + vV) | 

+ + V^) + -^Aa/S-y {SaS^S^ - SSaP/sP-y 



X exp < d x 



(3.5) 

Here we have introduced 18 scalar, pseudoscalar bosonic flelds aa, tTq, together with a vectorlike 
fleld V {a € {0, . . . ,8}) and, additionally, 18 auxiliary fields Sa,Pa necessary to deal with the 
six-fermion interactions induced by the 't Hooft term plus one auxiliary field V. In momentum 
space, the fermion determinant det ^ reads 

^{p,p') := {p\J\p') 

= {-f + m,)(27r)4 <5(4)(p -p')+e (^f^ A„ [a^{p - p') + i^^7:^{p - p')] ^3^^ 
-(p + p'\ f + f' , 

The major difference compared with the two-fiavor partition function (2.10) are the cubic ex- 
pressions in the auxiliary fields Sa,Pa- As a consequence an analytic evaluation of the functional 
integral over the auxiliary fields is not possible anymore. A suitable method to proceed at this 
point is the stationary phase approximation (SPA): choose the fields Sa,Pa so as to minimize 
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the integrand in the bosonized partition function (3.5). A necessary condition imposed on the 
fields is therefore: 

3H 

Ua + GSa + —^Aafi-f [>S's5'^' - P^P-y] = , 

3H (^-^^ 

TTq + GPa ^Aa/^jSijPj = , 

where Sa, Pa are now to be considered as (imphcit) functions of aa,T!'a- Performing, in addition, 
the integration over the vectorHke fields V the bosonized action can thus be written: 

= _ In det ^> ^ J d^x 

/( G H ^ C^-^) 

Prom here on we can apply the methods developed for the two-flavor case (Sects. 2.2-2.5) in order 
to reproduce the pseudoscalar meson spectrum within the nonlocal three-flavor NJL model. 

3.2 Mean-field approximation, gap equations and chiral condensates 

Starting from the action 5^°*^, Eq. (3.8), a power series expansion is performed around the 
expectation values of the fields (Tq, tTq, (a = 0, 1, . . . , 8) and v 

aa{x) = o-Q, + 5aa{x) , 

Tlaix) = 6T^a{x) , (3.9) 

v{x) = V + 5v{x) . 

Given the constraint imposed by charge conservation, one needs to consider only the fields 
(To,(T3,(78. (In the isospin limit studied later, one has the additional constraint that (73 also 
vanishes.) It is useful to introduce 

(7 = diag((7„, ad, cts) ■= o-qAq + craAa -I- asXs , (3.10) 

and analogously, S = diag(S'„, Sd, Ss) = SqXq + S3X3 + SgXs- Since {wa) = {Pa) = to leading 
order, the action in mean-field approximation reads 



cMF p J 4 




Trln [z-\p)(pH3y,s + M\p) 



+ ^-2] X] ((^iSi + ^SiSij+—SuSdSs 



(3.11) 



i&{u,d,s} 

where M{p) = diag (M„(p), Md(p), Ms(p)) with 

Miip) = Z{p)(mi + aiCip)) . (3.12) 

Here Z{p) is the wave-function renormalization factor (2.17), IsxS denotes the unity matrix in 
flavor space and F^^^ is the four-dimensional Euclidean volume. 

The mean-field equations (gap equations) are deduced once again by applying the principle 

of least action, = for cTj = ai {i E {u, d, s}), and — ^— = for v = v. The Si and Pj are 

both implicit functions of ai , determined through the SPA equations in mean-field approximation 
(compare Eq. (3.7)). One finds: 

- H — — 

CTj = —GSi '^^ijk^ijk^j^k (3.13a) 

^ «Ar f , Z{p)Mi{p) 

^--^ E / 7^ ^(.) (3.13c) 
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Finally, the chiral condensate (qq) is calculated in a way analogous to the two-flavor case 
using the definition (2.22). This leads to 

*> = -^~«/| 

now with q = u,d,s. Note that Mq{p) — )• niq for large p. As in Eq. (2.23), the subtraction of 
the free quark propagator makes sure that no perturbative artifacts are left in (qq) for niq 7^ 0. 



Z{p)Mq{p) 

p'^ + M^ip) 



p^ + 



(3.14) 



3.3 Masses and decay constants of the pseudoscalar meson nonet 
3.3.1 Second-order terms and meson masses 

We focus now on the masses of the (pseudoscalar) meson octet, plus the rj' meson with its 
large mass related to the anomalous breaking of the U(1)a symmetry. Consider once again 
second-order corrections to the mean-field action, extracted from a functional Taylor expansion. 



:(2) 



= \ I st^s '^^"^'^^ ^''^^''^ + 1 1 ^'"^ '^^'^^^^ ^""^^y^ 



5aa 5a 13 

?2 



where the second derivatives, 



S^Se 



S'S^ and 



are evaluated 



5aa(x)5<Tfj(y)'' S-Ka{x)SiT^{y)^ 5v(x)5v{y) Saa{x) 5v(y)- 

at the mean- field values o-„(.t) = CTq,, etc. Mixed second-order derivatives including one differ- 
entiation with respect to TTa{x) vanish because of tr{75} = trl'j^'y^} = tr{7^ 7,^75} = 0. 

A basis change with iVij = {Xa'^a)ij gives the standard representation of the pseudoscalar 
octet: 



^ V8 ^ VP 



TTij = (it) 




'78 I VO 



K+ \ 

2>?8 I M 
' V6^ V~3 



(3.15) 



Defining analogously a corresponding matrix a for the scalar boson octet, the fermion determi- 
nant, Eq. (3.6), can be written as 



^{P,p') = i-p + m,)(27r)^ 5(^)(p -p')+C (^^) V2 [a{p - p') + i75^(p - p')] 

'p + p'\ ^ + / 



(3.60 



2k 



-V{P-P') 



The resulting second-order contributions to the action are given by 
(2) _l f d^p 



+ 



1 



(2^)4 

d^p 

(27r)4 



with 



where 



nf^ip) = -8N, 



d\ 

(2^ 



^ijM^P^ ^'^vi.P) ScTkei-p) + Qij^keiP) ^T^ij{p) ^T^u{-p) 
[Qv{p) Sv{p) Sv{-p) + Qij,v{p) S(Tij{p) Sv{-p)] , 

q+ ■ q- T Mi{q+)M^{q-) 



e\q)Z{q+)Z{q 



[{q^Y + Mf{q+)\ [{q-f+MKq-)] 



(3.16) 



(3.17) 



(3.18) 
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5=*= = g lb |, and (r='=) ^ is given as the solution of the system 



(3.19) 



For completeness we have also determined the (inverse) propagators associated with the vector- 
type field v: 



1 4/Vc 



i&{u,d,s} 



q+-q- - Mi{q+)Mi{q~) + 


' (q^f (q-f - (q^ ■ q-f 






\q+f + Mf{q+)_ 


\q-f + Mf{q-) 





8Nc f d^g 



K J (27r)4 



C{q)J^iq)Z{q+)Z{q- 



q-[q-Mi{q+) + q+Mj{q-)] 



{q+f + Mf{q+) (q-f + MUq-) 



The meson masses can now be determined by writing the second-order term of the action, 
Eq. (3.16), in the physical basis as 



:(2) 



1 f d'^p 



{a.(/) [7r°(p) 7r°(-p) + 27r+(p) 7r-{-p)] + 



P 2 J (27r)4 

+ GKip^) [2K\p) K\-p) + 2K+{p) K-{-p)] + 

+ QssiP^) Vsi'P) Vsi-p) + Gooip'^) Vo{p) Voi-p) + 2^08 (P^) Voip) Vsi-p)} , 

where the inverse propagators Gp arc defined in Eqs. (3.17)-(3.19). Considering only the isospin- 
symmetric case, m„ = rrid with (73 = 0, one has 



H 



-1 



Gn{p')=(G + -Ss] +n-,{p') 
GKip') = (G + ^s^) +nzs{p' 



Gssip^) 
Gooip"^) 
Gosip"^) 



+ i7-Jp2^ + 2i7-(p2 



1 r 6G-4HSu-2HSs , ^- . 2 
3 [2^2 - GHSs - H^Sl ""^-^ 
1 r 6G + AHSu- HSs , 2^ rj-( 2x" 

3 [2G2 - GHSs - H^Sl + ""^^ ^ + \ 

H{Ss - Su) < JJ- ( 2\ _ TT- / 2n 

3 2G2 - GHS, ~ H^S?, ""^^ ' ''^^ ' 



(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 



The masses of the pseudoscalar octet mesons are given by the poles of their propagators or, 
equivalently, 

Gp{-ml) = 0, for P G {tt, K, 77} . (3.25) 

Finally, the physical 77 and r)' mesons are identified after a basis change and introducing the 
mixing angle 9 = 9{p^): 

77 = 778 COS Or, - 7/0 sin Or, 

ri = 778 sin 0^/ -|- 770 cos 0^/ , 

where 0,, = 9{—m'^),6ri' = 9{—mf^,). With the inverse 77 and 17' propagators Gr] and Gr]', respec- 
tively, instead of ^00 j ^88, ^08 we obtain for the mixing angle 

26^08(^2) 



tan 26" (p^) 



Gooip^) - Gssip^ 



(3.27) 
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and: 



Evidently, the mass splitting m^, —m^ is given by twice the square root appearing in Eq. (3.29). 
3.3.2 Renormalization and decay constants 

As in the two-flavor case, renormalized fields^, <^(p) = Z<^^^^(^(p), are introduced so that the 
quadratic part of the Lagrangian can be written in the standard form 

Given the pole structure of the propagators Qp^{p^), one arrives at an explicit expression for 
the renormalization constants: 



" dp2 



, for P G {tt, K, r]} . (3.30) 



The pseudoscalar meson decay constants in the three-flavor case are defined as 

{o\Jl^io)mp)) = -^fafsp, ^ {o\Ji^m^^ip)) = iuzi/\^, (3.31) 

where Jxa.^^) ~ '4'{^)l^l5^'^i^) denotes the axial-vector current. As described in Sect. 2.5 
we have to gauge the nonlocal action in Eq. (3.1) connecting the fields by a Wilson line and 
introducing a set of axial gauge fields .4^ (a G {0 . . . ,8}). The only modifications compared 

to the two-fiavor calculation are the replacements r ■ ^ — > XaAa, ^af^a ^md r ■ tt — > AoTTq 
in the gauged fermion determinant (2.31). Following the lines of the two-fiavor calculation, we 
find for the matrix elements (3.31):^*^ 

, ^ S'^ In det 



A=0 

t=o 



+ 



2iAre(AL^Af + AfAi^)x 

— fir 1 ^ 1 n Q:iZ(q+)Z(q-)^[i{q-) 



xTvj [\aq,'^c(q-'-a) ^('^-^f " ^^^^ 

\J, '''' dq^ 2) [{qtf+Mf{qt)\[{q~^f + M]{q^)\ J 

+ 2iN,{>^iX=; + Xp^x 



^ 2 y [(g+)VM2(g+)][(g-)VM|(g^)]J 

(3.32) 



^Here ip{p) stands generically for any of the fields 7ra(p), . . . 

^°As in the two-flavor case, the integration over a can be performed analytically [25] 
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ruu = rud 


rrig 


G 


H 


K 


2.9 MeV 


69.3 MeV 


1.63 fm^ 


-0.77 fm^ 


5.28 GeV 



Table 3: Current quark masses (rriu = md,ms) and coupling-strength parameters G, H and k. 



with 

p _ p (3.33) 

The decay constants are derived from the expressions (3.32) and their definitions, Eq. (3.31), by 
contraction with p^^: 

^-1/2 

/a/3 = ^PMJAa(0)MP))^^ ' (3-34) 

evaluated at the corresponding mass p'^ = —fn^- Owing to the properties of the Gcll-Mann 
matrices and assuming isospin symmetry (i. e., = m^) one has^^ fa^ = Sa^fir for a € {1, 2, 3} 
and faj3 = Sap fx for a G {4, 5, 6, 7}. 

3.4 Pcirameters and results 

The enlarged set of parameters in the three-flavor case now includes, in addition, the strange 
current quark mass and the coupling strength H of the U(l)A-breaking (nonlocal) 't Hooft- 
Kobayashi-Maskawa interaction. The profiles of the nonlocality distributions, C(p^) and J-'(p^), 
are left unchanged, assuming that these distributions are governed by (flavor-independent) gluon 
dynamics, so flavor-symmetry breaking effects can be ignored in these quantities. With the 
parameters collected in Table 3, the results for chiral condensates, dynamical quark masses 
at = and properties of the pscudoscalar meson nonet including octct-singlct mixing are 
summarized in Table 4. The input current quark masses, m^^d — 3 MeV and ~ 70 MeV, are 
once again compatible with MS masses at a renormalization scale /it ~ 3 GeV and correspondingly 
large magnitudes of quark condensates at that scale. The resulting pscudoscalar meson-octet 
masses are generally within less than 2 % of their measured values, while the pion decay constant 
comes out just slightly smaller than its value fi^^ ~ 86McV in the chiral limit, about 10% 
below its empirical value. The result for the ri-rj' mixing angle, = —29.9°, agrees well with 
the empirical value from Ref. [26], 9 = -29.0°. 

Fig. 3 shows the dynamically generated quasiparticlc masses of the u- and s-quarks, M„(p) = 
Z{p)('mu + d'uC{p)) and Ms{p) = Z {p)(ms + ag C{p)) , where the flavor-independent distributions 
C{p) and J^{p) enter. Note the interplay of the chiral condensates Eq. (3.14), the scalar fields in 
Eqs. (3.13) and the wave-function renormalization factor Z{p) in determining those dynamical 
quark masses. 



^^This follows from the fact, that the summands in Eq. (3.32) can be written as 

+ y^X^:)Aij = 2 Re (aL^ ) . (3.35) 

Assuming m„ = one gets the stated result. Note, in particular, that for the two-flavor case, Aij is independent 
of i,j, hence leading the anti-commutators of the Pauli matrices, as given in Eq. (2.34). 

^^Note the different definitions of the ri-ri' mixing angle in this work and in Ref. [26]. The cited number 
6 = —29.0° has, however, already been translated to the definition, Eq. (3.26), of the mixing angle used in the 
present work. 
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V 




MM 


-0.290 


-0.309 


-1.953 


0.369 


0.583 





niK 








fK 






0.142 


0.503 


0.553 


0.957 


0.084 


0.100 


2.1° 


-29.9° 



Table 4: Output for condensates, dynamical quark masses, pseudoscalar meson masses, decay 
constants (all given in GeV) and 77-77' mixing angles resulting from the parameters listed in 
Table 3. 




p [GeV] 

Figure 3: Dynamical quasiparticle masses of the u- and the s-quark resulting from the three- 
flavor gap equations (3.12). 

4 Thermodynamics: Polyakov loop and nonlocal PNJL model 

4.1 Prepcirations 

The nonlocal NJL approach provides a basic understanding for the mechanisms governing spon- 
taneous chiral symmetry breaking and the dynamical generation of massive quasiparticles from 
almost massless current quarks, in close contact with QCD itself. But the thermodynamics 
based on the NJL model goes wrong: colored quasiparticles would exist in regions of tempera- 
ture T and baryon chemical potential /^b where they should be suppressed by confinement. This 
is the point where the Polyakov loop enters. It acts as an order parameter for the confinement- 
deconfinement transition in pure-gauge QCD.^^ The coupling of the dynamical quarks to the 
Polyakov loop promotes the NJL scheme to the PNJL model [6-13, 16]. Once this coupling is 
active, it prevents colored objects (quarks, diquarks and color-octet qq pairs) from propagating 
in the hadronic low-temperature phase, while color-singlet (mesonic) quark-antiquark pairs and 
three-quark compounds are permitted (although the formation of localized baryons is not yet 

^^Alternative possibilities to describe the interplay between chiral symmetry breaking and confinement within 
the framework of NJL models exist in the literature (e.g., external chromomagnetic fields mimicking the gluon 
condensate [27]). 
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realized at this stage). 

In the PNJL context the Polyakov loop $ is expressed in terms of homogeneous (constant) 
background fields, and A^, representing the Euclidean temporal gauge- field components 
associated with the color-diagonal generators of the SU(3) gauge group: 



^ = ^ trc 

Nr. 



e^p(^-^{AlXs + AlXs)^ 



(4.1) 



Here A?c = 3 is the number of colors, trc indicates the trace over color space only, /3 := 1/T 
is the inverse temperature, and A3,A8 are the Cartan-algebra elements of the color SU(3) Lie 
group. One can derive an effective potential, the Polyakov-Fukushima potential ^(<I>, T), 
that describes the first-order confinement-deconfinement transition in the pure-gauge case (here 
we adopt the version proposed in Refs. [9, 10]): 



rp4 



--62(T)$*$-F64(r)ln l-6$*$-h4(^$*^-F$3^ -3($*$)2 . (4.2) 

This effective potential satisfies the condition of being symmetric under the center Z{3) of the 
gauge group: $ — t- exp[27rin/3] with n = 1,2,3. The temperature-dependent coefficients 
b2,A{T) are given in Ref. [9], parametrized such that the first-order deconfinement transition in 
pure-gauge lattice QCD [28] is accurately reproduced together with its transition temperature 
To = 270 MeV. Although Z//($,$*,T) is parametrized in terms of the longitudinal components 
A^, A\ of the Cartan algebra, it implicitly includes effects of the non-Cartan elements of color 
SU(3), as well as effects of the transversal components of the gluon fields A^^, as pointed out in 
Refs. [13,29,30]. 

The PNJL model introduces the coupling of the Polyakov loop to the quarks in the standard 
gaugc-covariant derivative and adds the effective potential U to the mean- field action (2.16). 
The thermodynamics is described using the well-known Matsubara formalism. The Euclidean 
time component p4 of the quark four-momentum is replaced by 

, q As . Q Ah 

P4 -cOn + + Ai— + Al— , 

with Matsubara frequencies LOn = (2ra -|- l)7rT, n G Z and quark chemical potential fx. For 
antiquarks, the corresponding replacement is 

PA -oJn - iM - Y + 
Four-momentum integrations are replaced as 



(2vr)4 ■ (27r) 



4.2 Two-flavor nonlocal PNJL model 

The PNJL thermodynamic potential with Nf = 2 quark flavors becomes^^ 



d'V 

2 (2^ 

nez ^ ' 



Tr In 



^(ia;„,p) 



+ 



-I- 



2G 



+ Z^($,$*,T). 



(4.3) 



Here 5 ^ is the inverse quark propagator expressed in Nambu-Gor'kov space, 





(Z Mia;„7o-7-p-M-i(^4-Fi/i)7o 
S '(iw„,p) = ^ 





(Z*)-i(ia;„7o-7-p-M*+i(A4+i/x)7o 



(4.4) 



^*Note an extra factor 1/2 because of the doubling ol the degrees of freedom in Nambu-Gor'kov space. This 
notation is useful for a separate treatment of particle and antiparticles. 
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where the momentum-dependent mass matrix M is diagonal in color space, 



M = diag, (M (uj- ,p),M{uj+,p),M{uIp)), 



(4.5a) 



and Z is the wave-function renormalization matrix, also diagonal in color space, 

Z = diag,(Z(a;-,p), Z{u+,p), Z{u;lp)) . 



(4.5b) 



In Eqs. (4.5a) und (4.5b), M{p4,p) is the dynamically generated mass, Eq. (2.18), and Z[p4,p) is 
the wave-function renormalization factor, Eq. (2.17), already encountered in the zero-temperature 
treatment. The generalized Matsubara frequencies appearing in M{p4,p) and Z{p4,p) are: 



± 



i^l ± Al/2 - Al/{2V3) , 



i/i + Al/Vs , 



(4.5c) 



with iOn = (2n + l)7r T,n G Z the standard fcrmionic Matsubara frequencies. 

In the limit of vanishing quark chemical potential, = 0, the thcrmodynamical potential 
(4.3) for N{ = 2 can be written in the simpler form 



n 



-AT 



E E 

ie{o,±} nez 



d^p 
(2^ 



In 



a;r+p' + M2(^;,p; 



+ 



2G 



+ W($,$*,T), (4.6) 



where the generalized Matsubara frequencies are again given in Eq. (4.5c). The finite- temperature 
generalization of the gap equations (2.21) are obtained by the requirement of a minimal grand- 
canonical potential O, leading to the necessary conditions 



dv 



dn on 



dAl 



dAl 



0. 



(4.7) 



According to Refs. [9, 10] one has $ = <I>* in mean-field approximation and, consequently, 
= 0. In the following we first discuss mean-field solutions at /i = 0. Eq. (4.7) determines the 
expectation values of the scalar and vector fields, a{T) and v{T), together with the Polyakov 
loop $(r). The mean fields a and v enter in the calculation of the thermal condensate, {qq) (T) = 
—dfl/dmo, via the dynamical quark mass M{u}l^,p) and the renormalization factor Z{ujl^,p), 
when taking the derivative of the pressure P = —i} with respect to the quark mass mo. 

A key result is shown in Fig. 4. As in our previous work [11], the entanglement of the chiral 
and deconfinement transitions is a characteristic result of the PNJL approach, a feature that 
persists in the nonlocal version including wave-function renormalization effects. The crossover 
transition temperature in this example is Tc = 205 MeV, to be compared with the critical 
temperature, Tq = 270 MeV, of the first-order deconfinement transition in pure-gauge lattice 
QCD that has been used to fix the temperature scale in the Polyakov-loop effective potential 
(4.2). In Fig. 4 no change has been made in this potential when introducing the coupling to 
dynamical quarks. In Ref. [31] it has been argued that the transition temperature, Tq, for the 
deconfinement transition depends on the number of dynamical quark fiavors: ro(A^f = 2) = 
(208 ± 30) MeV instead of To{Nf = 0) = 270 MeV mentioned earlier. Using this shift in To as 
input in the Polyakov-loop effective potential U, the chiral transition temperature is lowered to 
Tc ~ 170 MeV shown in Fig. 5. 

Compared to the nonlocal PNJL model without wave-function renormalization, there is now 
the condensate v that becomes temperature dependent. The result is shown in Fig. 6. It is 
evident, that v stays constant up to a temperature that is about 1.5rc. Above this temperature 
the absolute value of the condensate v decreases. The reason for this can easily be understood: 
the wave-function renormalization factor approaches one at high momenta (see the lattice data 
in the right picture of Fig. 1); a similar effect must appear when the temperature is increased. 



^^Recall that we are considering the isospin-symmetric limit m„ = md = mo- 
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TIT, 

Figure 4: Solid curves: calculated temperature dependence of the chiral condensate {qq) (left) 
and of the Polyakov loop $ (right) normalized to the transition temperature Tc = 205 MeV as 
obtained in the nonlocal PNJL model considered here. The dashed lines show the chiral conden- 
sate for the pure fermionic case and the Polyakov loop for the pure gluonic case, respectively. 



1.0 




T [GeV] 

Figure 5: The Tq dependence of the chiral condensate. The chiral transition tempera- 
ture decreases from 205 MeV to about 170 MeV when reducing To(A^f = 0) = 270 MeV to 
ro(iVf = 2) = 208 MeV. 



Finally, the pressure P = — O is computed after subtracting a divergent vacuum (T = 0) 
term. The result is shown in Fig. 7. This figure displays, in addition, the separate contributions 
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a. 




Figure 6: Temperature dependence of v. The 
decrease of (the absolute value of) v at high 
temperatures stems from the fact that the 
wave-function renormalization function Z{p) 
approaches one at high energy scales. 



Figure 7: Temperature dependence of the nor- 
malized pressure P/T'^ (at zero-density) (solid 
line). The dashed and dashed-dotted lines 
show the contributions of the Polyakov loop 
and the quarks to the pressure, respectively. 



from the Polyakov loop (i.e., the term AP = —lA in Eq. (4.6) featuring $(r) in the presence 
of quarks; dashed line) and from the quark quasiparticles (dashed-dotted line). Of course these 
separate contributions do not have a physical interpretation when taken individually, but their 
pronounced variations in the transition region give an impression of the strong coupling be- 
tween Polyakov loop and quarks around Tc- The total pressure is then a result of nontrivial 
cancellations between both parts as shown by the smooth solid curve in Fig. 7. 

4.3 Three-flavor nonlocal PNJL model 
4.3.1 Thermodynamical potential 

The thermodynamics of three quark flavors introduces the strange quark, with its heavier mass 
rus well separated from the masses m„ ^ of the light quarks. It also introduces the axial U(l) 
anomaly as a mechanism that generates additional nonlinear dynamics in the resulting gap 
equations [11, 12]. 

By analogy with the two-flavor case, the stationary phase approximation gives the grand- 
canonical mean- field thermodynamical potential^^ in terms of the fields (3.13): 



if / r \ H 1-2 ^^'^^ 

I fe{u,d,s} ^ ^ ) 



where the scalar mean fields now carry quark flavor indices. The inverse Nambu-Gor'kov quark- 
propagator matrix is 



S ^{\Un,p] 



'z ^(|iw„7o-7 •p-M-i(A4-|-i/x)7oj 

(i*)-i (ia;„7o-7 • p-M*+i(A4+iA)7o) 

(4.9) 



^®The factor \ results again from the doubling of the degrees of freedom in Nambu-Gor'kov space. 
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where the momentum-dependent dynamical mass matrix M is again diagonal in color and flavor 
space, 



/diagc (M„ (a;„_„ , p ) , M„ (a;+„ , p ) , M„ (0;° „ , p ) ) 



M = 



diag,(Mrf(a;-„,p), Mrf(a;+„,p-'), Mdls*] ,„,p)) 

diag,(M,(a;-„,p), M,(c^+„,p), M,(a;0„,p))y 

(4.10a) 

with Mi{p),i G {u, d, s} from Eq. (3.12). The matrix of the wave-function renormalization 
factors is 

Z = diag,(Z(a;-,p), Z{u+,p), Z{luIp)) , (4.10b) 

with 

i^% = ^n-ifif±Al/2-Al/{2V3), ul^ = Un-ifif + Al/V3. (4.10c) 

In contrast to the two-flavor model, different quark chemical potentials ft = diag^d^u, ^^di l^s) 
have now been introduced for the three quark flavors. The trace in Eq. (4.8) can be further 
simplified leading to 

2 



^ = -2^ E EE 

fe{u,d,s}i=0,±neZ 

1 { 



d'^p 
(2^ 



In 



' 1 / -J 

( fe{u,d,s} 



G 



H 



(4.11) 



E ( ^f^f + 9 ^/^/ ) + -^SuSdSs \ + ^+ ^(^, ^*,T). 



This is the thermodynamic potential of the nonlocal three-flavor PNJL model in mean-field 
approximation, including quark wave-function renormalization. The auxiliary scalar fields Sf are 
determined by the SPA conditions, Eq. (3.13a). Note that the Polyakov-loop effective potential 
U is taken to be universal, independent of quark flavors. 



4.3.2 Gap equations in mean-field approximation 

Given the thermodynamic potential (4.11), the fields (T„ = ad,(Ts, v and ^4,^4 are determined 
by requiring J7 to be stationary, as in the two-flavor case. The necessary conditions are given by 
the gap equations 

together with the stationary phase approximation equations (3.13a). As long as we limit our- 
selves to the zero-density [ji = 0) case we have ^1 = (cf. Ref. [9]). 

A first instructive result, displayed in Fig. 8, concerns the temperature dependence of the 
Polyakov loop <I) in the presence of u-, d- and s-quarks. The dashed curve in this figure shows 
the pure-gauge first-order deconfinement transition at the critical temperature Tq = 270 MeV. 
This phase transition is controlled by the effective potential (4.2), constructed to satisfy Z(3) 
symmetry and parametrized so as to reproduce the energy density of purc-gluc lattice-QCD 
thermodynamics. The coupling to quarks induces two basic effects. First, Z{2>) symmetry 
is explicitly broken by the presence of quarks and the first-order phase transition turns into 
a continuous crossover. Secondly, the transition-temperature range is lowered to a broad band 
between 150 and 200 MeV (see the solid curve in Fig. 8). This feature is also consistently observed 
by different lattice-QCD groups in their simulations with 2-t- 1 quark fiavors and the light quarks 
taken at (almost) physical masses [32,33]. Note that the effective potential U has been left 
unchanged when arriving at the (solid) crossover curve in Fig. 8, i.e., the parametrization of 
the coefficients in lA still maintains the pure-gauge temperature Tq = 270 MeV. The shift of the 
transition temperature from Tq to Tc below 200 MeV emerges entirely through the coupling of 
the dynamical quarks to the Polyakov loop in this calculation. 
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T [GeV] 

Figure 8: Results for the Polyakov loop in the presence of two light and one heavier quark 
compared to lattice-QCD computations (open circles: Ref. [33], filled squares: Ref. [32]). The 
solid line shows our PNJL model calculation with Tq = 270 MeV as input in the effective 
potential U. The dashed curve corresponds to the pure gauge case, compared to lattice data 
from Ref. [34]. 



At this point our findings appear to be at variance with considerations in Ref. [31] suggesting 
that the feedback induced by the dynamical quarks should change the Polyakov-loop effective 
potential such that ro(A^f = 0) = 270 MeV is shifted downward to ro(iVf = 2 + 1) = (187 ± 
30) MeV. Such a shift in U would run the crossover profile for the Polyakov loop (solid curve in 
Fig. 8) down to much lower temperatures in the present approach, out of scale for comparison 
with the A^f = 2 + 1 lattice-QCD data. On the other hand, a more sophisticated treatment of 
the effective potential lA using renormalization- group methods [35, 36] might help clarifying the 
situation. 

Fig. 9 shows the temperature dependence of the chiral condensates [uu) = {dd) and {ss) 
resulting from the present calculation. The softening of the strange quark condensate clearly 
reflects the stronger explicit chiral symmetry breaking by the s-quark mass m^. One notes that in 
these condensates, the influence of the quark wave-function renormalization factor Z{p) through 
the gap equations (3.13) and (3.14) is only marginal in comparison with previous calculations 
[11, 12] that used Z{p) = l. 

Recent (2-|-l)-flavor lattice-QCD computations of the chiral condensate have been performed 
using physical light quark masses by the Budapest-Marseille- Wuppertal (BMW) group [33, 37] 
and by the "hotQCD" collaboration [32, 38] (the latter employing a highly improved staggered 
quark (HISQ) action). These lattice data suggest a relatively smooth chiral crossover with tran- 
sition temperature in the range ~ (160-170) MeV. Our mean-field nonlocal PNJL calculation, 
using a Polyakov-loop effective potential U with To{Nf = 0) = 270 MeV as before, is not able 
to reproduce these updated lattice results (see dashed curve in Fig. 10). Running the Tq scale 
in U down to To(A^f = 2 + 1) = 187 MeV and adding pionic contributions from chiral pertur- 
bation theory [39] give the qualitatively correct pattern for {qq){T) (see solid curve in Fig. 10). 
However, as already mentioned, the Polyakov loop would now be systematically shifted to lower 
temperatures, out of scale with the corresponding lattice data. In addition, such small values of 
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To enhance the inference of the Polyakov loop such that the chiral transition tends to become 
first order, in contradiction with the lattice results. Again, improved Polyakov-loop effective po- 
tentials with systematic implementation of the backreaction from light quarks [35, 36] are under 
discussion in this context. 




^"8.0 0.1 0.2 0.3 0.4 

T [GeV] 

Figure 9: Results of nonlocal PNJL calculations of the chiral condensates including quark wave- 
function renormalization. Solid curve: temperature dependence of the chiral condensate {uu) = 
{dd). The strange-quark condensate {ss) is shown as the dashed curve. 



4.3.3 Finite densities and phase diagram for 2 + 1 flavors 

Consider now the phase diagram emerging from the nonlocal PNJL model at nonzero quark 
chemical potential. The starting point is the (2 + l)-flavor thermodynamical potential (4.11) 
treated in mean-field approximation. For calculational simplicity we choose the symmetric 
combination = /i„ = /Lt^ = /ig. The case ^Ug = together with = jid has been studied earlier 
in Ref. [12]. The flavor-symmetric choice of the quark chemical potential shifts the /U-scale just 
slightly upward as compared to the //^ = case. 

The resulting {T^jj) phase diagram. Fig. 11, is the generic one for PNJL type models at 
mean-field level: a chiral crossover transition ending in a critical point followed by a first-order 
transition line down to ~ 0.3 GcV at T = 0. The behavior of the Polyakov loop, detached from 
the chiral first-order-transition boundary has led to an interpretation in terms of an intermediate 
"quarkyonic" phase [40] . These features are not changed by the additional effects of quark wave- 
function renormalization in the nonlocal PNJL model. 

However, recent two- flavor investigations beyond mean field [35, 36], including both the quan- 
tum fluctuations in the quark quasiparticle sector and the backreaction from the matter fluc- 
tuations to the Polyakov-loop effective potential, suggest that the coincidence of chiral and 
deconfinement crossover transitions continues all the way down to a critical end point that now 
appears at low temperature, Tcep ~ 25 MeV, at /i ~ 0.3 GeV. Similar tendencies have been 
found in Ref. [41]. That is then the terrain of nuclear physics, however, where the PNJL ap- 
proach is not supposed to work and different methods (e.g., chiral effective field theory with 
baryons) must be applied [42]. 
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Figure 10: Solid curve: temperature dependence of the chiral up-quark condensate using 
To = 187 MeV and adding pionic contributions as obtained in chiral perturbation theory (e. g. 
Ref. [39]). Dashed curve: mean-field result with Tq = 270 MeV (same as solid curve in Fig. 9). 
Lattice-QCD data from Refs. [32,33,37,38]. 



Along the chiral crossover transition at small chemical potential fx, it is convenient to intro- 
duce a transition band, Tc{fi), specified by the region in which the chiral condensate has dropped 
to about half of its vacuum value: 0.4 < {ij^ip) / {Tpip) q < 0.6. This is also the domain in which 
the chiral susceptibility develops a pronounced maximum. In this region an expansion 

_P_.,_..(^)Vo((^)') ,4.13, 

is meaningful. From Fig. 11 one finds 

K = 0.11 ±0.02 (4.14) 

where the error refiects the uncertainty in locating the transition line in the crossover region. This 
should be compared to recent values found in lattice QCD: ii'iat — 0.06 from Ref. [43] (analysis 
of the scaling properties of the chiral condensate and its susceptibilities), and Xiat — 0.08 from 
Ref. [44] (Taylor-series expansion), both with an uncertainty of about 10%. 



5 Summary and outlook 

In this work the two- and three-flavor PNJL models have been further extended with inclu- 
sion of quark wave-function renormalization effects. This brings the PNJL approach in close 
contact with Dyson-Schwinger calculations and lattice-QCD results for the Landau-gauge quark 
quasiparticle propagator, with the wave-function renormalization factor Z{p) of the quark quasi- 
particle propagator introduced in addition to the momentum-dependent dynamical fermion mass 
M{p). 

In the following we summarize the main results of the present studies and point out improve- 
ments compared to previous work: 



25 




yu [GeV] 

Figure 11: Phase diagram for the (2 + l)-flavor nonlocal PNJL model at mean-field level. 
The orange band shows the confinement-deconfinement crossover transition as described by 
the Polyakov loop in the range 0.1 < <I> < 0.3. The dashed black line corresponds to the chi- 
ral crossover (blue band: 0.3 < ('0V')/('0V')o < 0.7). The solid black line indicates the chiral 
first-order transition. The temperature scale is set by = 0.2 GeV. 



(i) Including wave-function-renormalization effects requires a careful re-assessment of chi- 
ral low-energy theorems. Pseudoscalar meson masses and corresponding decay constants 
at zero temperature have been re-derived. The results clearly show that the formal- 
ism incorporates fundamental chiral relations such as the Gell-Mann-Oakes-Renner and 
Goldberger-Treiman relations. In the three-flavor case, the inclusion of the 't Hooft- 
Kobayashi-Maskawa interaction leads to the correct mass splitting between the r] and the 
r]' meson. 

(ii) The PNJL thermodynamics has now been developed with systematic inclusion of the quark 
quasiparticle renormalization factor Z[p). The temperature dependence of the chiral con- 
densate and of the Polyakov loop has been calculated, indicating chiral and deconfinement 
crossover transitions. We have compared our results with recent lattice-QCD computa- 
tions. Finally, a quark chemical potential has been introduced that enables extensions to 
the finite-density region of the QCD phase diagram. 

(iii) The impact of the wave- function renormalization factor Z[p) compared to previous calcu- 
lations [11, 12] setting Z[p) = 1 is generally quite small over the whole relevant momentum 
range. This can be understood considering the gap equations at zero temperature: since 
Z{p) deviates significantly from unity only in the momentum range p ^ 1 GeV, its ef- 
fect does not contribute much to the relevant integrals because of its suppression by the 
integration measure. 

(iv) With inclusion of Z{p)^ the chiral and deconfinement crossover transitions tend to become 
smoother compared to our previous investigations. 

(v) The fiavor dependence, To(A'^f), of the deconfinement temperature scale is an important 
issue in lowering chiral-transition temperatures in accordance with the tendency recently 
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observed in lattice-QCD computations. We have investigated the impact of different values 
for To in the Polyakov-Ioop effective potential and found that such modifications do not 
relax the strong entanglement between chiral condensate and Polyakov loop, at least at 
mean-field level. 

(vi) Concerning the behavior of the phase diagram at finite quark (or baryon) chemical po- 
tentials jj, the backreaction of the matter fluctuations on the Polyakov-loop sector is an 
important effect, as pointed out in Refs. [35, 36], that calls for further studies. Finally, the 
challenge remains to construct a phase diagram at finite baryon densities in accordance 
with realistic constraints from nuclear physics. 
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